Infrared behavior of gluon and ghost propagators from asymmetric lattices 
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We present a numerical study of the lattice Landau gluon and ghost propagators in three-dimensional 
pure SU (2) gauge theory. Data have been obtained using asymmetric lattices ( V = 20 2 x 60, 40 2 x 60, 
8 2 x 64, 8 2 x 140, 12 2 x 140 and 16 2 x 140) for the lattice coupling /3 = 3.4, in the scaling region. 
We find that the gluon (respectively ghost) propagator is suppressed (respec. enhanced) at small 
momenta in the limit of large lattice volume V . By comparing these results with data obtained 
using symmetric lattices (V = 60 3 and 140 3 ), we find that both propagators suffer from systematic 
effects in the infrared region (p < 650MeV). In particular, the gluon (respec. ghost) propagator is 
less IR-suppressed (respec. enhanced) than in the symmetric case. We discuss possible implications 
of the use of asymmetric lattices. 

PACS numbers: 11.15.Ha, 12.38.Aw, 12.38.Lg, 14.70.Dj 



I. INTRODUCTION 



Gluon and ghost propagators are powerful tools in 
the (non-perturbative) investigation of the infra-red (IR) 
limit of QCD and are thus of central importance for un- 
derstanding confinement (see e. g- D1J2D- In fact, ac- 
cording to the Gribov-Zwanziger IE and to 
the Kugo-Ojima ja, 0] confinement scenarios in Landau 
gauge, the ghost propagator must show a divergent be- 
havior — stronger than p~ 2 — for vanishing momen- 
tum p. This strong infrared divergence corresponds to 
a long-range interaction in real space, which may be 
related to quark confinement. At the same time, ac- 
cording to the former scenario, the gluon propagator 
must be suppressed and may go to zero in the IR limit 
HHHBQ, EHHd This would imply that the 
real-space gluon propagator is positive and negative in 
equal measure, i.e. reflection positivity is maximally vio- 
lated 0,0], indicating absence of gluons from the phys- 
ical spectrum (gluon confinement). These predictions are 
also valid for the case of pure SU{2) gauge theory and 
for three space-time dimensions. 

The nonperturbative evaluation of gluon and ghost 
propagators can be achieved by first principles methods 
such as Dyson-Schwinger equations (DSE's) 0,0 or the 
numerical simulation of lattice QCD. Landau-gauge stud- 
ies of DSE's (see e. g. ^IHjEI) have found an IR be- 



havior of the form D(p ) 



P 



[implying 



-D(O) = if k > 0.5] for the gluon propagator and of the 
form G(p 2 ) ~ £,- 2k - 2 = p-2a G -2 f or g^Qgt propa- 
gator, with the same exponent k (i.e. with cld = 2ac)- 
The available predictions for the IR exponent point to- 
wards k > 0.5 for pure SU(N C ) gauge theory in the 
four dimensional case. For the 3d case the exponents 
are « 0.4 and cld ~ 1.3. Note that in the d dimen- 
sional case ^| the relation between cld and aa is 
given by cld = 2ac + (4 — d)/2, implying for the quan- 
tity a(p 2 ) = (g 2 /4Tr)D(p 2 )G 2 (p 2 )p 6 the infrared behavior 
p2(a D -2a G ) — p4-d_ Thus, in the 4-dimensional case the 
running coupling cx(p 2 ) displays and infrared fixed point. 



Numerical studies of lattice gauge theories confirm the 
IR divergence of the Landau ghost propagator |Isl . lTgll20l 
I2l1 ] and have now also established that the Landau gluon 
propagator shows a turnover in the IR region, attaining 
a finite value for p ss 0. (A reliable extrapolation of D(0) 
to the infinite- volume limit is still lacking H^]-) More 
precisely, indications of a decreasing gluon propagator for 
small p have been obtained in the 4d SU(2) and SU (3) 
cases for the strong-coupling regime Il9l.l2a.l24l . in the 3d 
SU{2) case (also in the scaling region) |25l 1261 127| . in the 
3d SU{2) adjoint Higgs model [H and in the 4d SU(2) 
case at finite temperature [28|. The actual turning of 
the gluon propagator has been seen clearly for the equal- 
time three-dimensional transverse gluon propagator in 4d 
SU(2) Coulomb gauge [llllljor the 3d SU{2) Landau 
case using very large lattices (23 (of 140 lattice sites per 
direction) and — recently — in the Ad S U(3) Landau 
case with the use of asymmetric lattices 113, H3, 0] • 

Thus, the two nonperturbative approaches above seem 
to support the Gribov-Zwanziger and the Kugo-Ojima 
confinement scenarios in Landau gauge. However, the 
agreement between the two methods is still at the qual- 
itative level, at least when considering the gluon and 
the ghost propagators. Moreover, recent lattice stud- 
ies |35l |36( seem to indicate a null IR limit for a{p 2 ), 
instead of a finite nonzero value. Of course, in this 
case, special care must be taken to eliminate finite- 
size effects, especially when the IR region is considered 
IH El El El IH OS. This issue has also been stud- 
ied analitically using DSE's |3^| , showing that torus and 
continuum solutions are qualitatively different. This sug- 
gests a nontrivial relation between studies on compact 
and on noncompact manifolds. Finally, the nonrenor- 
malizability of the ghost-gluon vertex — proven at the 
perturbative level plj \ confirmed on the lattice [36l l4l| 
(for p > 200 MeV) and used in DSE studies to simplify 
the coupled set of equations — has been recently criti- 
cized in Ref. |35| . Thus, clear quantitative understanding 
of the two confinement scenarios is still an open problem. 

In this work we extend the study presented in [2^ for 
the 3d SU(2) case, by including results for the ghost 
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propagators from very large lattices and by conside ring 
also asymmetric lattices, as done in |32> H3, E3 13 
\a!\ . Our aim is to verify possible systematic effects re- 
lated to the use of asymmetric lattices (as suggested in 
|36|). by comparing the results to the ones obtained for 
symmetric lattices. In particular, we focus on the de- 
termination of the IR critical exponents ao and ag and 
on their extrapolation to the infinite- volume limit. The 
study of gluon and ghost propagators in three space-time 
dimensions is computationally much simpler than in the 
four-dimensional case and it may help to get a better 
understanding of the Ad case Note that the 3d 

case is also of interest in finite-temperature QCD (see for 
example |4fiLl47j V 

In the following section we describe our simulations 
and present our results. We end with our conclusions. 



II. SIMULATIONS AND RESULTS 



s 




Simulations have been done using the standard Wilson 
action for SU (2) lattice gauge theory in three dimensions 
with periodic boundary conditions. We consider /? = 3.4 
and several lattice volumes, up to V — 140 3 . In Ta- 
ble^ we report, for each lattice volume V, the parame- 
ters used for the simulations. All our runs start with a 
random gauge configuration. We used the hybrid overre- 
laxed algorithm (HOR) 48, 49] for thermalization and the 
stochastic overrelaxation algorithm or the so-called Cor- 
nell method for gauge fixing thermalized configurations 
to the lattice Landau gauge |50l lip ), |52l . The numeri- 
cal code is parallelized using MPI [53. For the random 
number generator we use a double-precision implemen- 
tation of RANLUX (version 2.1) with luxury level set to 
2. Computations were performed on the PC clusters at 
the IFSC-USP. We ran part of the 60 3 lattices on three 
nodes and the 140 3 lattices on four or eight nodes. 

In order to set the physical scale we consider 3d SU (2) 
lattice results for the string tension and the input value 
\J~o = 0.44 GeV, which is a typical value for this quantity 
in the Ad SU(3) case. With H = c = 1, this implies 



TABLE I: Lattice volumes V, configurations, HOR sweeps 
used for thermalization and between two configurations (used 
for evaluating both propagators). 



V 


Configurations 


Thermalization 


Sweeps 


20 3 


500 


550 


50 


20 2 x 40 


800 


440 


40 


20 2 x 60 


600 


550 


50 


60 3 


400 


1100 


100 


8 2 x 64 


100 


550 


50 


8 2 x 140 


1000 


550 


50 


12 2 x 140 


600 


660 


60 


16 2 x 140 


600 


770 


70 


140 3 


101 


1650 


150 



FIG. 1: Plot of the gluon propagator D(p 2 ) as a function 
of p (both in lattice units) for lattice volumes V = 20 3 (x), 
20 2 x 40 (□), 20 2 x 60 (O) and V = 60 3 (Q). 



1 fm~ l = 0.4485^/er. For (3 = 3.4 we obtain the average 
plaquette value (W^i) = 0.672730(3). Then, considering 
the tadpole-improved coupling (3i = /3(Wi,i), we can 
evaluate the string tension y/b~ in lattice units using the 
fit reported in Eq. 2 and Table IV of Ref. 0. (The fit 
is valid for > 3.0, i.e. the coupling (3 considered here is 
well above the strong-coupling region.) This gives ^fo = 
0.506(5), implying a lattice spacing a = 0.227(2) fm, i.e. 
a -1 = 0.869(8) GeV. Also, we are able to consider here 
momenta as small as p m i n = 39.0(4) MeV and physical 
lattice sides of almost 32 fm. Finally, let us notice that, 
if we compare the data for the string tension (in lattice 
units) with the data reported in Ref. [55| for the SU (2) 
group in four dimensions (see also Ref. |21|). then our 
value of (3 corresponds to (3 w 2.21 in the Ad case. 

We study the lattice gluon and ghost propagators [re- 
spectively, D(p 2 ) and G(p 2 )] as a function of the mag- 
nitude of the lattice momentum p{k) (see Ref. [27|, H(j 
for definitions). We consider all vectors k = (k x , k y , kt) 
with only one component different from zero. For asym- 
metric lattices, the non-zero component has been taken 
along the elongated direction; for symmetric lattices, 
when evaluating the gluon propagator, we average over 
the three directions in order to increase statistics. Here, 
we do not check for possible effects of the breaking of ro- 
tational invariance Nevertheless, as we will see 
below, systematic effects due to the use of asymmetric 
lattices are evident mostly in the IR limit, where one ex- 
pects the breaking of the rotational symmetry to play a 
small effect (see also Ref. p2|h 

In Figures 2] [21 and [31 we plot the data obtained for the 
gluon propagator for different lattice volumes. (Errors 
represent one standard deviation.) One can clearly see 
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TABLE II: Fit of the gluon dressing function Z(p 2 ) = D(p 2 )p 2 
using the fitting function z (p 2 ) a ° (l + aip 2 + (i2p 4 ) and data in 
the interval p 2 < 0.85/a 2 . For each lattice volume V we report 
the infrared exponent an, the number of degrees of freedom 
(d.o.f.) and the x 2 /d.o.f . For the lattice volume V = 20 2 
there were not enough data points to do the fit. 

V a D x 2 /d.o.f. #d.o.f. 



20 2 x 40 


1.026 (7) 


0.4 


2 


20 2 x 60 


1.028 (5) 


0.4 


5 


60 3 


1.133 (9) 


1.1 


5 


8 2 x 64 


0.94 (4) 


1.2 


5 


8 2 x 140 


0.981 (5) 


1..3 


17 


12 2 x 140 


0.991 (5) 


1.2 


17 


16 2 x 140 


1.019 (6) 


1.2 


17 


140 3 


1.13 (1) 


1.2 


17 



FIG. 2: Plot of the gluon propagator D(p 2 ) as a function of 
p (both in lattice units) for lattice volumes V = 8 2 x 140 (x), 
12 2 x 140 (□), 16 2 x 140 (O) and V = 140 3 (Q). 
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FIG. 4: Plot of the gluon IR exponent an as a function of 
the inverse lattice volume 1/V (in lattice units) for lattice 
volumes: V = 20 2 x 40 and 20 2 x 60 (□), V = 8 2 x 64, 
8 2 x 140, 12 2 x 140 and 16 2 x 140 (O), V = 60 3 and 140 3 
(x). Note the logarithmic scale on the x axis. 



FIG. 3: Plot of the gluon propagator D(p 2 ) as a function of 
p (both in lattice units) for lattice volumes V — 8 2 x 140 ( x ) 
and 20 3 (□). 



that the propagator is IR-suppressed for sufficiently large 
lattice volumes. This behavior is evident both for the 
asymmetric and the symmetric lattices. However, the 
shape of the propagator is clearly different in the two 
cases for momenta p < 0.75/a, i.e. p < 0.65 GeV. Indeed, 
for asymmetric lattices the gluon propagator starts to 
decrease only at very small momenta (see Figures ^ an d 
|5J|. A similar behavior can also be observed in Fig. 2 of 
Ref. [13 and in Fig. 4 of Ref. 34]. Thus, considering 



the asymmetric lattices one would estimate a value M < 
0.25/a w 0.22 GeV as a turnover point in the IR, i.e. 
the momentum p = M for which the propagator reaches 
its peak. On the other hand, considering the two largest 
lattice volumes, i.e. V = 60 3 (see Fig. [IJ and V = 140 3 
(see Fig.0), the gluon propagator is clearly a decreasing 
function of p for p < 0.5 /a, corresponding to M < 0.435 
GeV. This is in agreement with Refs. |22t |27| where the 
value M = 350^5°° MeV is reported. We thus see a 
difference of almost a factor 2 between the momentum- 
turnover point in the symmetric and asymmetric cases. 

Let us also note that in Figures^and|21the gluon prop- 
agator at zero momentum D(0) is monotonically decreas- 
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FIG. 5: Plot of the ghost propagator G(p 2 ) as a function of 
p (both in lattice units) for lattice volumes V — 8 2 x 140 ( x ) 
and 20 3 (□). Note the logarithmic scale on the y axis. 
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FIG. 6: Plot of the ghost propagator G{p 2 ) as a function of p 
(both in lattice units) for lattice volumes V = 8 2 x 140 (x), 
12 2 x 140(D), 16 2 x 140(O) and V = 140 3 (CO- Note the 
logarithmic scale on the y axis. 



ing when considered as a function of the lattice volume 
V. On the other hand, if one considers the symmetric 
lattice V = 20 3 = 8000 and the slightly larger asymmet- 
ric lattice V — 8 2 x 140 = 8960, one finds (see Figure 
|3J) that D(0) is larger (by about 20%) in the asymmetric 
case. Thus, an extrapolation of D(0) to infinite volume 
using data from asymmetric lattices is most likely also 
affected by systematic effects. 



TABLE III: Fit of the ghost propagator G(p 2 ) using the fitting 
function z/((p 2 ) 1+aa (l+aip 2 +a,2p 4 )) and data in the interval 
p 2 < 2 /a 2 . For each lattice volume V we report the infrared 
exponent ac the number of degrees of freedom (d.o.f.) and 
the value of x 2 /d.o.f . 

V a G x 2 /d.o.f. #d.o.f. 



20 3 


0.272 (6) 


0.4 


1 


20 2 x 40 


0.18 (1) 


4.9 


6 


20 2 x 60 


0.154 (9) 


3.9 


11 


60 3 


0.21 (1) 


4.6 


11 


8 2 x 64 


0.053 (7) 


1.7 


12 


8 2 x 140 


0.025 (2) 


1.9 


31 


12 2 x 140 


0.060 (4) 


2.0 


31 


16 2 x 140 


0.093 (6) 


2.0 


31 


140 3 


0.174 (7) 


8.5 


31 



We also tried a fit to the data using the fitting functions 
considered in Ref. [34| (see results in Table HIT) . We see 
that the IR exponent a d usually increases with the lattice 
volume, i n ag reement with the results reported in Table 
2 of Ref. |2J|. However, it is clear from Figure 01 that 
the dependence of ajj on l/V is not simple. Indeed, old 
looks almost constant when considering V = 20 2 x 40 
and 20 2 x 60, but it gets a larger value when using the 
symmetric — and much larger — lattice volumes V = 
60 3 and 140 3 . A smaller value for oo is also found if 
one tries an extrapolation as a function of l/V for the 
results obtained using the three asymmetric lattices V — 
x 140. Again, it could be difficult to have systematic 
errors under control when trying an extrapolation using 
only asymmetric lattices. 

Strong systematic effects can also be observed in the 
ghost case. The propagator is less IR divergent when one 
considers asymmetric lattices (see Figures and EJ ■ As 
a consequence, the IR exponent a G is also smaller (see 
Table ITTTf for asymmetric lattices. This is the case either 
when comparing an asymmetric lattice to a symmetric 
one with almost equal lattice volume (i.e. V — 8 2 x 140 
and V = 20 3 ) or when considering an asymmetric lattice 
and a symmetric one with the same largest lattice side 
(e.g. V = 16 2 x 140 and V — 140 3 ). Moreover, the depen- 
dence on the lattice volume is very different in the two 
cases. Indeed, for symmetric lattices (e.g. V — 20 3 , 60 3 
and 140 3 ), the IR exponent decreases as the lattice vol- 
ume V increases (at fixed /?). This is in agreement with 
the result obtained in the Ad SU(2) case (see Table 4 and 
Figure 3 in Ref. [l^)- O n the other hand, a G increases 
as a function of V for the lattices V = -/V 2 x 140 with 
increasing N s . Thus, there is no simple relation between 
a G and V and an extrapolation in V using asymmetric 
lattices is probably of difficult interpretation. 

Using the above results we can evaluate the quantity 
a D - 2a G - 1/2, which is zero using DSE's [Hll6j. We 
find a D - 2aa ~ 1/2 = 0.21(2) for V = 60 3 and 0.28(2) 
for V = 140 3 . For the asymmetric lattices this value is 
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somewhat larger. However, let us note that these results 
depend on the choice of the fitting function. Moreover, 
the null value for arj — 2ac — 1/2 should be obtained 
only in the infinite-volume and in the continuum limit. 
We will analyze these limits in a future study. 

In this work we did not do a careful stu dy o f Gribov- 
copy effects [HEiE3,ElElElilElElE3, since we 
are interested in possible systematic effects due to the use 
of asymmetric lattices [6!|. Nevertheless, for the gluon 
propagator and V = 20 3 and 8 2 x 140, we compare results 
obtained using a standard gauge-fixing algorithm to data 
obtained with the so-called smearing method |4lLl68l|. (In 
this case we considered 2000 configurations for each lat- 
tice volume and each gauge- fixing method.) We found 
that for V = 20 3 the value of D(0) is 0.873(9) when us- 
ing the smearing method, while one gets 0.906(9) with a 
standard gauge-fixing method. The two results differ by 
several standard deviations and the propagator is smaller 
when using the smearing method, i.e. as one gets closer 
to the so-called fundamental modular region. On the 
contrary, for V = 8 2 x 140 we found D(0) = 1.183(15) 
with smearing and D(Q) = 1.152(15) without smearing. 
In this case the relative difference between the two cases 
is smaller and the result using smearing is larger than 
in the standard case. This seems to be a systematic ef- 
fect, which we observe for the first nine smaller momenta, 
when considering the asymmetric lattice. We believe that 
this issue deserves a more careful analysis. 

III. CONCLUSIONS 

We have compared data for gluon and ghost propaga- 
tors (in minimal Landau gauge) obtained using symmet- 



ric and asymmetric lattices. We find, for both propa- 
gators, clear evidences of systematic effects at relatively 
small momenta (p < 650MeV w 1.5^/a). In particu- 
lar, the gluon (respectively, ghost) propagator is less sup- 
pressed (respec. enhanced) in the IR limit when consid- 
ering asymmetric lattices than for the case of symmetric 
lattices. This implies that the estimates for the IR criti- 
cal exponents ac and od are systematically smaller in the 
asymmetric case compared to the symmetric one. Also, 
for the gluon propagator, the turnover point M is sig- 
nificantly smaller when considering asymmetric lattices 
than for the symmetric ones. Let us recall that M corre- 
sponds to the Gribov mass scale in a Gribov-likc propa- 
gator D{p 2 ) = p 2 /(p 4 + M 4 ). Finally, we have seen that 
the extrapolation to infinite volume of results obtained 
using asymmetric lattices is also most likely affected by 
systematic errors. We conclude that, even though using 
an asymmetric lattice does not modify the qualitative be- 
havior of the two propagators, one should be careful in 
extracting quantitative information from such studies. 

Our data have been obtained in the 3d SU(2) case. 
However, the behavior observed for the gluon propagator 
is very similar to what is obtained in the Ad SU(3) case 
pSlt l32l l33l |23 , where of course a study of this type is 
more complicated. 
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